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A computer-based technique has been developed to search bright 
BATSE gamma-ray bursts for spectral lines in a comprehensive man- 
ner. The first results of the search are discussed and an example line 
candidate shown. 



MOTIVATION 



^— I ■ Prior to the launch of CGRO, and during the early days of the mission, 
K>" ' the BATSE team expected that GRB lines would be common and easily idcn- 

. tified. Our expectations were based upon the observations of KONUS (|l|), 
^ ; HEAD A-4 (|l|) and Gmga (00). Our early efforts were based upon these 

i expectations: we searched bright bursts visually, examining consecutive spec- 
T— I I tra and also the spectrum of the entire burst (Eoh. Our expectations were not 

i met: no significant lines have been found (p^,|5|)T 

■ While the failure to find lines in BATSE GRBs does not yet imply a serious 
f~| . discrepancy between the Ginga and BATSE results (p^,p|), it is disappoint- 
On' ing. Previously, lines were interpreted as strongly supporting the theory that 
Q . GRBs originate nearby from a disk population. Results from BATSE on the 

' spatial distribution of GRBs demonstrate that only a minority of the bursts 
C/3 . can originate from a disk population ([I^ jfj,0) , so if the existence of lines is 

■ confirmed, their properties might indicate the presence of a disk subclass or 
J> . give clues to the physics of halo or cosmological sources. 

Is the failure to detect lines in the BATSE data because lines do not exist 
or is the failure due to some inadequacy in the visual search approach? To 
' answer these questions a more systematic search is needed. We have therefore 
implemented an automatic, computer-based comprehensive line search. The 
goal of the search is comprehensiveness — a brute force approach is taken so 
that no significant line will be missed. The purpose of the search is to identify 
line candidates — the search need not perfectly evaluate their significances, 
which will be done later under direct human control. The computer-based 
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search has several advantages in addition to comprehensiveness: subjectivity 
is ehminated, there is no variation in detection threshold due to human ex- 
haustion, there is no bias towards absorption or emission features, and the 
search will collect statistics on all trial lines. 

COMPREHENSIVE SEARCH METHOD 

The addition to the BATSE instrument of the eight Spectroscopy Detectors 
(SDs) was largely motivated by the KONUS results. Simulations (ppl|) and 
tests ( p^ of the performance of the BATSE SDs show that they should be 
capable of detecting KONUS or Ginga-\ike GRB lines. Each SD is a Nal(Tl) 
scintillator crystal, 12.7 cm diameter and 7.6 cm thick, viewed by a photo- 
multiplier tube of the same diameter. To limit the effort to a manageable 
level, the initial searches are restricted to lines below 100 keV, so only detec- 
tors in high-gain mode are useful. Typically, such detectors cover the energy 
range of about 10 keV to about 1400 keV, although analysis begins at about 
15 or 20 keV due to an electronic artifact (|l|). Work to model this artifact 
continues, so that in the future we will be able to extend the analysis to lower 
energies. Analysis of the data is based upon an instrument response model 
which accounts for the detectors, nearby spacecraft material, and scattering 
from the Earth's atmosphere ([2C|). When a burst occurs, high-time resolu- 
tion SHERB (Spectroscopy High Energy Resolution Burst) data are collected 
using a time-to-spill algorithm. 

Since we do not know a priori when or how long a line will exist, we search 
essentially all time scales available. We search each single SHERB record, 
every consecutive pair, triple, and group of 4, 5, 6, 8, 10, 14, 20, ... records 
and the sum of all the SHERB records. Clearly, many overlapping intervals 
are searched, and thus the searched intervals are not all independent. Addi- 
tionally, many intervals will have insufficient signal-to-noise ratio for a real 
line to be detectable — these intervals serve as controls. While not every pos- 
sible consecutive combination of records is formed, enough combinations are 
searched so that no line should be missed, although its significance might not 
be optimized. 

Similarly, we do not know a priori at what energies lines occur. Therefore, 
we first fit each spectrum with a continuum model, and then we perform 
continuum plus line fits using a closely spaced, fixed grid of trial centroids. 
Since the trial centroids are separated by one-third the detector resolution 
FWHM, no line candidate should be missed, although the exact centroid will 
not be found by the search. 

Candidates are identified by large 20) changes in x^- We have switched 
to using Ax^ instead of the F-test because it is more appropriate when the 
errors are known (^4|J|,^). After a candidate is identified, other time intervals 
will be tried and the centroid will be made a free parameter in the fits. 

To increase the robustness of the automatic nonlinear fits, we use a con- 
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FIG. 1. Results of the comprehen- 
sive hne search for GRB 920627 
(trigger 1676), SD 2. The fig- 
ure shows the change in result- 
ing from adding to the continuum 
model a narrow, additive Gaussian 
line. While SD 2 has useable data 
from 15-1240 keV, the search was 
made using the Comptonized con- 
tinuum model and the restricted en- 
ergy range 15-390 keV. The few 
points with Ax^ < are due to 
poor convergence. 



tinuum model with a small number of parameters, namely the Comptonized 
model, which is a power-law times an exponential cutoff. This is fit to a 
restricted energy range of « 400 keV. Our initial search is confined to lines 
which are narrow compared to the detector resolution and we use the simplest 
line model: an additive (or subtractive) Gaussian. So that there is continuum 
at energies below any line, the lowest energy trial centroid is two detector- 
resolution FWHM above the starting energy of the fits. Since very few of the 
fits begin below 15 keV, only rarely do we search for 20 keV lines. 

We show as an example of a burst without a significant line the data of 
SD 2 for GRB 920627. The 63 SHERB records of SD 2 are summed into 
693 spectra which are fit with trial lines at 20 different centroids from 27.3 
to 100 keV. The largest change in obtained by adding a line was 13.3 
(Fig. 1), which has a chance probability in a single fit of a few tenths of a 
percent. Considering the large number of spectra searched (albeit not all 
independent), this is insignificant. 

FIRST RESULTS OF THE COMPREHENSIVE SEARCH 

The search results for GRB 940703 are quite different from those for 
GRB 920627. Only one high gain detector, SD 5, viewed the burst at an 
angle less than 100°. The data from SD 5 comprises 79 SHERB records, from 
which 917 spectra were formed. The largest value of Ax^ identified in the 
search was 58.2 and many overlapping intervals had Ax^ values approach- 
ing this value (see Fig. 2). When more careful fits are made, using the Band 
spectral form (g) and the entire available energy range, the candidate is found 
to have Ax^ of "only" 23.4. The difference in the values of Ax^ is because 
for this very bright burst the Comptonized model is inadequate even over a 
restricted energy range. 

Most of the line significance appears to come from a shorter interval which 




20 40 60 80 100 
Trial Centroid (keV) 



4 



60 



FIG. 2. Results of the comprehen- 
sive Hne seareh for GRB 940703 
(trigger 3057), SD 5. The fig- 
ure shows the change in x'^ result- 
ing from adding to the continuum 
model a narrow, additive Gaussian 
line. While SD 5 has useable data 
from 20-1380 keV, the search was 
made using the Comptonized con- 
tinuum model and the restricted en- 
ergy range 20-397 keV. The angle of 
the detector normal with the burst 
was 51° and with the geocenter, 
79°. The few points with Ax^ < 
are due to poor convergence. 
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excludes periods of weak emission, for which Ax^ = 28.1. The spectrum 
of this interval is shown in Fig. 3 — the candidate is seen to be an emission 
feature with an equivalent width of 2.1 keV. The line is narrow — allowing the 
width to vary does not significantly reduce x^- For 3 additional parameters, 
the chance probability in any particular fit of such a large Ax^, if no line 
actually exists, is 3 x 10~^. Additionally, the residual plots show that adding 
the line to the model eliminates conspicuous nms of high, then low, residuals. 
Even considering the large number of spectra searched, the line is highly 
significant -the number of effectively independent spectra is well below 917 
and the line is significant in more than one spectrum. 

To date, 42 of the brightest GRBs have been searched by this technique. 
There was an average of 2.1 detectors per burst and 53 SHERB records per 
detector, for a total of 4729 SHERB records. From these records 51,651 
spectra were formed and 861,372 fits performed. Partial examination of the 
results has identified 8 candidates with Ax^ > 20. There are both absorption 
and emission candidates and their centroids range from 40 to 70 keV. 

More work remains to be done. In most cases several high-gain SDs ob- 
serve a gamma-ray burst. We can therefore usually test the reality of a line 
candidate by examining the data of the other detectors. If the lines are real, 
in all cases the data should be consistent and in some cases there should be 
confirmation: statistically significant detections in more than one detector. 
The ability to perform these tests is a major advantage of BATSE, but it 
is also time-consuming because there are many issues. Until we have com- 
pleted these tests, the BATSE team considers the features identified by the 
comprehensive search to be line candidates rather than detections. 
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FIG. 3. Data from SD 5, GRB 940703, for the interval 42.112-64.704 s. Top 
panels: Count rate data and fits, shown up to 100 keV. Bottom panels: Count 
rate residuals (data — model) in units of a for the entire energy range used in the 
fits. Left panels: A continuum-only fit, using the Band spectral form; = 212.0 
with 199 degrees- of- freedom. Right panels: Added to the fits is an narrow, additive 
Gaussian line; ~ 183.9 with 197 degrees-of-freedom. The "hump" at 30 keV in 
the data and continuum model is expected from the detector physics (^,|l8|). The 
contribution of atmospheric scattering to the counts in the line region is 
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